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ABSTRACT 

The latest Sloan Digital Sky Survey data reveal a prominent bifurcation in the distribution of debris 
of the Sagittarius dwarf spheroidal (Sgr) beginning at a right ascension of a w 190°. Two branches of 
the stream (A and B) persist at roughly the same heliocentric distance over at least 50° of arc. There 
is also evidence for a more distant structure (C) well behind the A branch. This paper provides the 
first explanation for the bifurcation. It is caused by the projection of the young leading (A) and old 
trailing (B) tidal arms of the Sgr, whilst the old leading arm (C) lies well behind A. This explanation 
is only possible if the halo is close to spherical, as the angular difference between the branches is a 
measure of the precession of the orbital plane. 

Subject headings: Galaxy: halo — Galaxy: structure — galaxies: dwarf — galaxies: individual Sgr 
dSph 



1. INTRODUCTION 

The disrupting Sagittarius dwarf spheroidal galaxy 
(Sgr) was discovered by Ibata, Gilmore & Irwin in 1994. 
It was soon realized that the S gr provided a po werful 
tool for the study of the Galaxy (Ib ata et alJll997tl . The 
nucleus of the Sgr has survived for many orbits around 
the Galaxy, whilst its tidal tails h ave now been detected 
over a full 360° on the Sky (see e.g.lTotten fc Irwinlll998t 
IMaiewski et al.l 120031 IBelokurov et alJ 12006(1 . The dis- 
rupted fragments of the Sgr diffuse in the Galactic po- 
tential. As the pericenter of the Sgr's orbit is ~ 16 kpc, 
whilst its apocenter is ~ 60 kpc, the debris provides a 
strong constraint on the Galaxy's halo. 

The morphology of the Sgr stream is known in de- 
tail in t he Galactic southern hemisphere thanks to 
2MA SS (IMaiewski et all 120031 12001 ISkrutskie et, all 
HHI. INewberg et all (120021 1200^ " made early detec- 
tions of Sgr tidal deb ris in data from the Sloan Digital 
Sky Survey (SPSS, seelHogg et al.l20 0lllStoughton et all 
[2002 ] ISmith et al) 120021: IPier et alJ I2003t llvezic et all 
12005 iGunn et al J 12006ft . Then. IBelokurov et alJ l)2006ft 
used a color cut to pick out the upper main sequence 
and turn-off stars in SPSS Data Release 5 (DR5, 
lAdelman-McCarthv et al.l2006ll belonging to the stream. 
In the Galactic northern hemisphere, they found a promi- 
nent bifurcation or branching in the stream, beginning at 
a right ascension a s» 190° (see the upper right panel of 
Figure The lower and upper declination branches of 
the stream, labelled A and B, can be traced until right 
ascensions of at least a as 140°. Using the location of the 
subgiant branch as an estimator, the A and B branches 
are reckoned to be at similar distances. There is also ev- 
idence in the data of a fainter, still more distant stream 
(C) directly behind the A branch. 
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Belokurov et al.'s (2006) dataset is important for two 
reasons. First, it traces the Sgr stream around the North 
Galactic Cap, the very spot at which obl ate and prolat e 
dark halos give different predictions (e.g., Helmi 2004a). 
Second, the debris in 2MASS is dynamically younger 
than that found in SPSS, so the SPSS data should give 
stronger constraints, as the stars have had longer to move 
in the Galactic potential. 

Early explorations of the evolution of the Sgr sug- 
gest ed that the Galactic halo may be close to spheri- 
cal dlbata et aUl2001bUMaiewski et all 1200311 . For ex- 
ample. IJohnston et al J l)2005jl showed that the precession 
apparent in Sgr debris in the 2MASS d ataset s trongly 
favored mildly oblate halos. However, Helmi (2004a) 
pointed out that many of the earlier datasets are re- 
stricted to stars that have only recently been torn off 
the Sgr and s o have n ot diffu sed in the Galactic poten- 
tial. In fact, iHelmil l|2004btl argued that the velocity 
measurements of 2MASS M giants in the leading arm 
favor strongly prolate halos an d this was subsequently 
confirmed bv lLaw et, alJ l)2005ll . Hence, present studies 
of the disruption of the Sgr have reached an impasse, 
with different datasets pointing to dramatically different 
flattenings. 

Here, our aim is to show how the newly discovered bi- 
furcation arises. Using numerical simulations, we argue 
that the complex morphology of the Sgr stream uncov- 
ered by Belokurov et al. (2006) can only be reproduced 
if the Galaxy's halo is close to spherical. Although our 
models do not resolve the contradiction between the pre- 
cession rate and the velocities in the leading arm, they 
do provide a new and powerful argument in favor of an 
almost spherical Galaxy halo. 

2. THE SIGNIFICANCE OF THE BIFURCATION 

The results of a typical simulation of the tidal disrup- 
tion of the Sgr in a nearly spherical potential are shown 
in Fig.^ We will give the details of the simulation set-up 
shortly, but at the moment our aim is to gain a quali- 
tative understanding of why the bifurcation occurs. The 
particles in Fig.Q] ar e color-coded according to when they 
were torn off the Sgr. In the direction of the SPSS PR5 
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Fig. 1. — Left: Simulation showing the tails of the Sgr dSph. Particles are color-coded according to when they were lost (gold: < 4 Gyr 
ago, red: between 4 and 5.7 Gyr, green: between 5.7 and 7.4 Gyr and blue > 7.4 Gyr ago). The yellow (orange) curves show the past 
(future) behavior of the Sgr's orbit over 2 Gyr. The orbital period is 0.7 Gyr. The positions of the Galactic centre (GC), the Sun, and 
Sagittarius (Sag) are marked. The green lines show the right ascension range 110° < a < 220°, which corresponds to the SDSS data 
analyzed by Belokurov et al. (2006). The 4 streams are marked A (young leading arm), B (old trailing arm), C (old leading arm) and D 
(young trailing arm). The circle gives the distance cut-off at 20 kpc. Upper Right: The SDSS data from Belokurov et al. (2006), with 
stars color-coded according to magnitude. Middle and Lower Right: Scatter plots in right ascension and declination of the tidal debris. 
Only particles within (beyond) a heliocentric distance of 20 kpc are plotted in the middle (lower) panels. The black squares show the field 
locations of Belokurov et al. (2006). In the middle panel, streams A and B are clearly visible. The upper arm is the old trailing material, 
while the lower arm is the young leading material. In the lower panel, the old leading material is in the lower and the young trailing 
material in the upper branch. [The simulation uses a Miyamoto-Nagai disk and logarithmic halo with = 1.05, together with the set d 
of proper motions. The mass of Sgr is 10 s Mq] 



TABLE 1 

Sets of proper motions used in the simulations. 



Label 


)l a COS <5 


He 


frad 


Remarks 




(in masyr - 


1 ) (in masyr - 


1 ) (in krns -1 ) 




a 


-2.65 


-0.88 


137 


HST measurement (Ibata et al 2001b) 


b 


-2.8 


-1.4 


137 


Schmidt plates measurement (Irwin et al. 1996) 


c 


-2.9 


-1.5 


137 


< 1<t variation of Schmidt plates 


d 


-3.02 


-1.49 


137 


Simulation fit from Law et al. (2005) 


e 


-3.05 


1.28 


137 


5cr variation of HST values 



data - namely, the opening angle denned by the green 
lines - there are 4 distinct streams of material. They 
are the young leading arm (labelled A), the old trailing 
arm (B), the old leading arm (C) and the young trailing 
arm (D). Here, old and young indicate when the stars 
were torn off. 3 out of the 4 streams are identified in the 
SDSS DR5 dataset. Stars belonging to the D stream are 
more difficult to detect as they occur in DR5 primarily 
in the range 180° < a < 220°, where they are not easy 
to untangle from the other streams. The simulation data 
are separated according to heliocentric distance and then 
projected onto the sky as viewed from the Sun, as shown 



in the lower right panels of Fig. ^ The young leading 
arm provides branch A and the old trailing arm branch 
B of the bifurcated stream of Belokurov et al. (2006). 
These two narrow branches are at similar heliocentric 
distances, as required to match the data. The material 
in these branches is about two revolutions apart in or- 
bital phase. For the material beyond 20 kpc, the older 
leading material is in the lower declination branch, while 
the younger trailing material is in the upper. The old 
leading arm (C) provides the more distant and fainter 
stream detected by Belokurov et al. (2006) behind the 
A branch. The young trailing arm (D) lies behind the B 
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Fig. 2. — The projected density of the young leading and old trailing tidal debris of the Sgr in a sequence of simulations with different halo 
Battenings, proper motions and Galaxy models. The key at the top of each column has the following significance: ML means Miyamoto- 
Nagai disk and logarithmic halo models, DB means the Dchnen & Binney models, whilst a, b, c, d or e refers to the choice of proper 
motions. The panels in a given column differ only in the flattening of the halo, with q^ for the ML models and q for the DB models recorded 
in the upper left corner. The mass of the Sgr is 10 8 Mq. Note that the particles in the simulation represent both stars and dark matter. 



branch and so is not the structure seen by Belokurov et 
al. (2006). 

The Sun lies roughly in the orbital plane of the Sgr. 
If the potential were exactly spherical, the debris of the 
Sgr would lie in a single plane and no bifurcation would 
exist. Any asphericity (whether intrinsic to the halo 
or produced by the bulge and disk) causes the orbital 
plane to precess and therefore the planes of the 4 arms 
to be slightly different. The positional difference between 
branches A and B is a direct measure of the precession 
over two orbital revolutions and hence the asphericity of 
the potential. The facts that (i) branches A and B are so 
close in projection and (ii) branch C lies behind branch 
A suggest that the precession is small, and that the po- 
tential is close to spherical. If the halo is too oblate or 
prolate, then debris is scattered over a wide range of loca- 
tions and does not lie in thin, almost overlapping streams 
on the sky. To back up this qualitative argument, let us 
now describe a suite of simulations developed to mea- 
sure the properties of the bifurcation as a function of 
halo flattening, Sgr mass and proper motion. 

3. SIMULATIONS 
3.1. Set-up 

The present position of the Sgr dSph is (a, 5) — 
(283.7°, —30.5°), while its heliocentric distance is 25 ± 2 
kpc and radial velocity is 137 kms -1 ijlbata et al.lll997j) . 
Listed in Table ^ are two mea surements of the p roper 
motion of the Sgr, the first from llr win et al.l l)1996j) using 



TABLE 2 

The Strength and Range of the Bifurcation in Nearly 
Spherical Haloes 



Logarithmic Halo a 



Q<t> 




a 


b 


c 


d 


e 


0.92 


1.0 




1.7 


1.9 


2.0 


1.5 








170° 


180° 


190° 


160° 


0.95 


1.05 




2.7 


2.3 


2.4 


1.7 








190° 


190° 


170° 


150° 



Dehnen & Binney Models' 3 



Q<t> 


1 


a 


b 


c 


d 


c 


0.95 


0.95 




1.1 


1.1 












160° 


150° 






0.97 


1.0 




1.4 


1.4 












130° 


150° 







a All Galactic models with q = 0.8, 0.9, 0.95, 1.0, 1.05, 1.11, 1.25 
and 1.5 were investigated. If there is no bifurcation, or if the lower 
branch of the bifurcation bends back to negative declinations, the 
model is discarded. For models in which there is a bifurcation, 
the strength (d), and the onset a a are given for each set of proper 
motions. 

b All Galactic models with q = 0.8, 0.9, 0.95, 1.0, 1.05, 1.11, 1.25 
and 1.5 were investigated but only those with bifurcations are re- 
ported. 



Schmidt plates and the second from llbata et alJ ((2001b) 
using HST data. Dines cu et al.'s (2005) re cent measure- 
ment agrees with that of llrwin et alJ l)1996|) within the er- 
rors. Given a choice of proper motions, we integrate back 
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Fig. 3. — Plots of the heliocentric distance and velocity versus 
right ascension for the same model as in Figure^ Again, particles 
are color-coded according to when they were lost (gold: < 4 Gyr 
ago, red: between 4 and 5.7 Gyr, green: between 5.7 and 7.4 Gyr 
and blue > 7.4 Gyr ago). The datapoints give the heliocentric dis- 
tances to the streams, as derived from fitting the subgiant branch 
described in Belokurov et al. (2006). Note that the distances to 
the A and B streams are too small for right ascensions a < 190°. 



in time for 10 Gyr adopting a potential for the Galaxy. 
At the final position, we insert a Plummer sphere con- 
taining 10 6 particles with a scalelength of 350 pc. We 
investigate models with a total mass of between 10 s and 
10 9 M©. The particles are integrated fo rward using the 
parti cle-mesh N-Body code Stjperbox (Fellha uer et al.l 
2000) until the position today is reached again. 

While the present position of the Sgr is unchanged in 
all our simulations, the range of proper motions recorded 
in Table Q] is investigated. For the Galactic potential, we 
use one of two possibilities. In the first (denoted by ML), 
the halo is represented by a logarithmic potential of the 
form 

$haioM-|ln(i? 2 + zV+^ 2 ), (1) 

with vo — 186 kms -1 and d = 12 kpc (where R and 
z are cylindrical coordinates). The parameter q$ is the 
axis ratio of the equipotentials. It controls whether the 
halo is spherical (g$ = 1), oblate (q$ < 1) or prolate 
(g$ > 1). In general, <?$ is of course not the same as 
the axis ratio in the density q, which varies with radius 
for the logarithmic potential (see e.g., Evans 1993). The 
disc is represented by a Miyamoto-Nagai potential: 



$disc(i?,z) = 



GM d 



R 2 + (b + Vz 2 + c 2 Y 



(2) 



with M d = 10 11 M , b = 6.5 kpc and c = 0.26 kpc 
Finally, the bulge is modelled as a Hernquist potential 

GM h 



$bulge(r) 



(3) 



using Mb = 3.4 x 10 10 M and a = 0.7 kpc. The su- 
perposition of these components gives quite a good rep- 
resentation of the Milky Way. The circular speed at the 
solar radius is ~ 220 kms -1 . The major advantage is 
the analytical accessibility of all quantities (forces, den- 
sities, and so on). Hence, this model has been very widely 
used - in particula r, in many previous investigations of 
a ( e 

ILawetalJia. 

In the second (deno ted by DB), we use the Galactic 
potential suggested bv lDehnen fc Binnevl l)1998j) . It con- 
sists of three disc components, namely the ISM, the thin 
and the thick disc, each of the form 



the Sgr stream ( e.g. lHehiill20f)4albl: ITohnston ^eTafll2005l: 



p disc (i?, z) = - — exp 
2z d 



Rm 
R 



R 

Rd 



(4) 



With i? m = 0, equation Q describes a standard double 
exponential disc with scale-length R d , scale- height z d and 
central surface-density Ed- For the stellar disks R m is 
set to zero, while for the ISM disc, we allow for a central 
depression by setting R m = 4 kpc. Furthermore, the halo 
and the bulge are represented by two spherical density 
distributions of the form 



ps(R,z)=p — 



, — 7 / \ 7-/3 

771 \ / 771 x 

1 + — 



'"(> 



exp 



m , 

--?-). (5) 



where m 2 = R 2 + z 2 q~ 2 and q is the axis ratio in the 
density. We ch oose the parameter s accor ding to the best- 
fit model 4 in iDehnen fc Binnevl l)1998j) . This provides 
a better representation of the Galaxy, but at somewhat 
greater computational cost. 

3.2. Results 

Snapshots of the distribution of tidal debris around 
the North Galactic Cap for some typical simulations are 
shown in Figure [3 We quickly see that most of the 
models do not look at all like the data. Only haloes 
close to spherical provide bifurcated streams. The simu- 
lated streams in moderately and strongly oblate or pro- 
late haloes do not bifurcate. 

To proceed further, we need to develop an objec- 
tive criterion for identifying the bifurcation. We use 
a Marquand-Levenberg routine to fit a single Gaussian 
and two Gaussians to the declination distribution in 
the young leading and old trailing tidal debris. Mod- 
els for which a single Gaussian is everywhere preferred 
(as judged by the A% 2 ) are unacceptable, as they do not 
show two identifiable streams. If two Gaussians are a 
better fit than a single, then the ratio d of the distance 
between the two peaks to the sum of the dispersions of 
each peak is computed as a function of right ascension. 
We refer to the mean value of this parameter, taken over 
all right ascensions 110° < a < 220°, as the strength of 
the bifurcation, (d). The onset of the bifurcation ccq is 
taken to be the right ascension when d = 1.5. Applying 
this algorithm to Belokurov et al.'s (2006) dataset, the 
bifurcation has strength (d) « 1.7 and begins at a right 
ascension ao ~ 190°. For a range of simulations, the 
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same quantities are recorded in Table [21 Both Galaxy 
models contain a flattened disk and bulge, so the model 
parameters {q^ and q) are not a reliable guide to the over- 
all flattening of the potential. Rather, we give in Tabled 
the axis ratio of the equipotentials at the mean of the 
pericentric and apocentric distances of the Sgr's orbit. 

There are a number of interesting conclusions from the 
Table. First, very few models actually give bifurcations 
at all. Of the 80 models tested, only 10 give bifurcated 
streams. A bifurcation occurs if the axis ratio of the 
potential lies in the range 0.92 < < 0.97. The 
best overall match to the data is given by simulations us- 
ing the Miyamoto-Nagai disk and logarithmic halo with 
q,jj = 1, together with sets c or d of proper motions. 
They reproduce the strength and the location of the bi- 
furcation reasonably well. Second, if we use the proper 
motion s measured by HST (set a, derived bv llbata et alJ 
we do not obtain bifurcated streams, whatever 
the halo flattening. The precession is controlled not just 
by the flattening of the potential, but also by the eccen- 
tricity of the orbit, and hence the proper motions. Sets 
b and c of proper motions, with values close to those 
derived from Schmidt p lates ijl rwin et aT1ll996|) . provide 
much better fits. Third, IHelmil <l2004l¥ has claimed that 
strongly prolate haloes with ijk 1.65 give the best fit to 
the 2MASS data on the right ascension, declination, he- 
liocentric distance and radial velocity of the tidal debris. 
Such strongly prolate models do not match the bifur- 
cated stream in the SDSS dataset. In fact, the claims of 
prolateness rely heavily on their choice of Galactic poten- 
tial (Miyamoto-Nagai disk and logarithmic halo) and are 
not reproduced with the Dehnen & Binney models. The 



Miyamoto-Nagai disk declines like a power-law, rather 
than an exponential, and so it is reasonable to inter- 
pret the finding of prolateness or stretching of the halo 
as compensation for deficiencies in the disk model. Simi- 
larly, although some of the nominally prolate halo models 
(<70 > 1) in Table [3 provide matches to the bifurcation 
data, the equipotentials are mildly oblate {Q^ < 1) at 
the radii probed by Sgr's orbit. 

Nearly spherical models can reproduce the projected 
density of the Sgr stream around the Northern Galactic 
Cap. However, as is traditional in this area, our simu- 
lations do not fit all the data! For example, there is a 
mismatch of ~ 20° in the right ascension of the beginning 
of the C stream (c.f., the upper and lower right panels of 
Figure Figure shows the heliocentric distances and 
velocities of stream A and B for the simulation of Fig- 
ure ^ together with the heliocentric distances derived 
by Belokurov et al. (2006) from subgiant branch fitting. 
The heliocentric distances of the simulated streams are 
within the observational error bars over the range of right 
ascensions a > 190°, but they are too small for a < 190°. 
In c ommon with other s imulations in oblate haloes (see 
e.g.. Uohnston et all2003l . the radial velocities of 2MASS 
M giants in the leading arm are also not matched. Pos- 
sible causes of these discrepancies are discussed shortly. 

The existence of the bifurcation also constrains the 
mass of the Sgr. For example, Figure 01 shows a sim- 
ulation identical to that of Figure ^ but with the initial 
mass of the Sgr increased to 5 x 10 8 M Q . The more mas- 
sive the Sgr, the greater its internal velocity dispersion. 
This has two consequences - the tidal arms are broader 
and they diffuse away from the Sgr's orbital path more 
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220 200 180 150 140 120 
RA [deg] 

Fig. 5. — A sequence of simulations differing only in the mass 
of the Sgr dSph at the beginning of our 10 Gyr simulations. This 
is marked in the top left-hand corner and varies from 2.5 X 1O 8 M0 
(top panel) to 10 9 Mq (bottom panel). If there is a bifurca- 
tion, then strength (d) is also recorded. The Galactic model is 
a Miyamoto-Nagai disk and logarithmic halo with = 1.05, while 
proper motion set d is used. 



TABLE 3 

Properties of the Models of the Sgr dSph. 

Initial Mass Final Mass Final Velocity Final Crossing 

(in 1O 8 M ) (in 10 8 M Q ) Dispersion (in kms" 1 ) Time (in Myr) 



1. 


0.5 


19. 


59 


2.5 


1.9 


30. 


38 


5. 


3.4 


36. 


46 


7.5 


5.7 


14. 


37 


10 


8.0 


50. 


32 



quickly. The overall effect is that debris is scattered over 
a wider range of locations. As the subpanels on the right 
hand side show, the bifurcation persists but is less dra- 
matic than in the data, whilst the A and B streams are 
no longer as collimated as in the data. The run of he- 
liocentric distances with right ascensions, however, is a 
somewhat better match. Figure shows a sequence of 
simulations with Sgr masses between 1O 8 M0 and 1O 9 M , 
while the Galactic potential is kept fixed as a Miyamoto- 



Nagai disk and logarithmic halo with q$ = 1.05. The 
bifurcation blurs with increasing mass. Quantitively, the 
strength of the bifurcation (d) falls from 2.4 when the 
Sgr's mass is 1O 8 M to 1.9 at 2.5 x 1O 8 M and to 1.2 
at 5 x 10 8 M Q . Once the Sgr mass rises much above 
5 x 1O 8 M , there is no visible bifurcation. High mass 
models are disfavored until it has been demonstrated that 
the tidal streams of stars can remain as highly collimated 
as in the data. In Table |3| we give the properties of the 
Sgr remnant at the end of our simulations. Of course, our 
simulations show both dark matter particles and stars, 
whereas the observable today is the luminous matter left 
in the Sgr dSph. The mass loss of the Sgr is mainly 
affected by its orbit and hence the choice of proper mo- 
tion together with the choice of Galactic potential. Since 
these are uncertain, we also give in Table [3] the internal 
crossing times at the virial radius and the three dimen- 
sional velocity dispersion, which governs the broadening 
of the tails. 

Our overall picture gives three predictions. First, 
a second wrap of branch D may be detectable in the 
2MA SS data, closer than that already reported by 
IMaiewski et al.l l|2003|) . Secondly, the dynamically older 
B stream should have a larger velocity dispersion than 
the younger A stream. The mean velocities of branches 
A and B probably differ by ~ 15 kms -1 , though this may 
be hard to measure as it may be less than the stream's 
internal dispersions. Thirdly, if our current ideas on the 
star f ormation history of the Sgr are correct fe.g.. lGrebel 
2000), then there may be a difference in the stellar pop- 
ulations of the streams. Stream A may contain evidence 
for a younger population that is not present in stream B. 

4. CONCLUSIONS 

The Sgr stream, as seen by the Sloan Digital Sky Sur- 
vey (Belokurov et al. 2006), is composed of two branches 
(A and B) at about the same heliocentric distance visibly 
diverging at right ascensions a ss 190° to give a bifurca- 
tion, together with a third stream (C) aligned with the 
A branch, but well behind it. This complex and intri- 
cate morphology throws down an enormous challenge to 
modellers. 

Here, we have given a physical picture of how this 
structure may arise. The bifurcation is caused by a pro- 
jection of the young leading (A) and the old trailing (B) 
tidal arms of the Sgr, while the old leading arm (C) lies 
well behind A. The bifurcation between A and B, and 
the positioning of C behind A, can only be reproduced 
in simulations of the disruption of the Sgr if the halo 
is nearly spherical. Simulations in either moderately or 
strongly oblate or prolate haloes fail these tests by a wide 
margin. In particular, the bifurcation only exists if the 
axis ratio of the potential Q</> a t the radii sampled by 
the Sgr's orbit lies in the range 0.92 < Q < 0.97. The 
physical explanation of this is easy to provide. The ma- 
terial in the A and B branches is about two revolutions 
apart in orbital phase. The angular difference on the 
sky between the A and B branches is therefore a direct 
measure of the precession of the orbital plane of the Sgr 
over two revolutions. As the angular difference is small, 
so the precession of the orbital plane is small, and so the 
potential muss to be close to spherical. If the potential 
is moderately prolate or oblate, debris is scattered over 
a much wider range of locations. The path from obser- 



The Halo Shape 



7 



vation to theoretical conclusion is surprisingly direct and 
independent of detailed modelling. 

The bifurcation also provides a strong constraint on the 
mass of the Sgr and its debris. If this is much larger than 
5 x 10 8 M Q , then the tidal streams are too diffuse to give 
a clear bifurcation. This also is easy to understand. As 
the internal velocity dispersion increases with progenitor 
mass, so the streams become broader and diffuse more 
quickly in the Galactic potential. The A and B branches 
do not then have the highly collimated appearance seen 
in the data. 

Our simulations - like all other simulations of the dis- 
ruption of the Sgr in the literature - do not agree with all 
the data. In particular, the detailed distances to the A, B 
and C streams given in Belokurov et al. (2006) are not re- 
produced over the full range of right ascension. Although 
this is a defect, the limitations of the commonly-used 
methodology for Sgr disruption simulations also need to 
be acknowledged. The underlying assumption is that the 
Galactic potential is static and unevolving over up to 
10 Gyr timescales. This is clearly incorrect - the Milky 
Way is believed to h ave accreted 30% of its mass over the 
last 5 Gyr (see e.g.. Ivan den Boscbll2002t fNeistein et al.l 
2006). Simulations typically show that the last major 
merger take place at about a redshift z = 1, roughly 
8 Gyrs ago l|Navarro. Frenk fc White! IT995h . Although 
most of the mass will have accreted in the outer parts of 
the Galaxy, the Sgr's orbit extends out to ~ 60 kpc and 
will surely been affected by this rearrangement. The for- 
mation of the Galactic bar has been dated to between 5 
and 8 Gyr ago (see e.g.. lSevensterlfl999albf) . and so bar- 
driven evolution of the inner Galaxy will also have caused 
substantial changes. The effects of time evolution are of 
much greater importance for the SDSS dataset than for 
the 2MASS dataset, which is restricted to dynamically 
younger material. 

The strength of the argument presented in this pa- 
per is that it relies on the gross morphological features 
of the Sgr stream. To reproduce the detailed positions 
and velocities of stars in the A, B and C branches may 
well require a clearer understand of Galactic evolution. 
However, the existence of a bifurcation in nearly spher- 
ical potentials is a robust result. To challenge the main 
conclusion of this paper requires the devising of an alter- 
native explanation of the existence of two streams that 
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are closely matched in distance over a ~ 50° arc. In this 
respect, our argument compares favourably with other 
methods of determination of halo shape using meth- 
ods such as the flaring of the neutral gas layer (e.g. 
Qllin g~ Merrifi gidl l2QQfl l) or the s tellar kinematics of 
halo stars fe.g.. lvan der Marellll991j) . These are afflicted 
by systematic uncertainties regarding the contribution of 
the cosmic ray pressure or the orientation of the stellar 
velocity ellipsoid, for example. In contrast, the bifurca- 
tion in the Sgr stream is a clean, simple and direct test. 
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